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ABSTRACT 

We study the survival of gas planets around stars with masses in the range 
1-5 Mq, as these stars evolve off the Main Sequence. We show that planets 
with masses smaller than one Jupiter mass do not survive the Planetary Nebula 
phase if located initially at orbital distances smaller than (3-5) AU. Planets more 
massive than two Jupiter masses around low mass (1 Mq on the Main Sequence) 
stars survive the Planetary Nebula stage down to orbital distances of ~3 AU. As 
the star evolves through the Planetary Nebula phase, an evaporation outflow will 
be established at the planet's surface. Evaporating planets may be detected using 
spectroscopic observations. Planets around white dwarfs with masses y^wD ^ 
0.7 M0 are generally expected to be found at orbital radii r > 15 AU. If planets 
are found at smaller orbital radii around massive white dwarfs, they had to form 
as the result of the merger of two white dwarfs. 

Subject headings: Stars:AGB and post-AGBs — Planetary Systems — white dwarfs 



1. INTRODUCTION 

In recent years, the quest for Jupiter-like giant planets has been extended to a com- 
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Friedrich et al.l 120051 ). The underlying assumption has been that planets can survive the 
parent star's evolution. Once this is accepted, the fact that white dwarfs are 10'^ to 10^ 
times fainter than their main sequence prog enitors opens up the possibility to observe plan- 
ets thr ough direct imaging in the infrare d (IBurleigh et al.ll2002l ) where the planet emission 
peaks (IBurrows et al.lll997l : Ilgnacell200ll ). 



In general, it has been assumed that planets survive to the white dwarf stage if they 
manage to stay in a large enough orbit to avoid engulfment by the star when the latter 
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increases its radius as it ascends the Red Giant Branch (RGB) and Asymptotic Giant Branch 
(AGB). Therefore, studies of the planet's fate as the star evolves off the main sequence have 



generally been restricted to determining the planet's orbit (e.g. iLivio fc Sokei 



1994: Sackmann et al. 1993: Rasioetal 



19961: 



Rvbicki fc Denisll200ll : iBurleigh et al 



1984 



Soker 



20021 ) 



or the orbit's stability (e.g. IPuncan fc Lissauei 19981 ). However, before the star reaches 



the white dwarf stage, and immediately after the AGB evolution, the star evolves into 
the planetary nebula (PN) phase, where the stellar temperature can reach 300 000 K with 
luminosities of 10^ Lq, while emitting powerful winds. 

In this paper we explore the effects of an evolving AGB and post-AGB star on an 
orbiting planet. In §2 we examine the planet survival during the AGB phase and the orbital 
changes due to the AGB mass- loss. In §3 we estimate the thermal conditions at the planet's 
surface as the star increases its temperature during the PN phase. In §4 we estimate the 
range of parameters under which an outflow from the planet will ensue, and we estimate the 
planet's evaporation rate due to thermal heating. In §5 we qualitatively discuss the planet's 
reaction to mass-loss and other processes that might influence the planet's evolution, and in 
§6 we discuss the planet's overall survival. Our conclusions follow. 



2. PLANET SURVIVAL DURING THE AGB PHASE 



Single white dwarf progenitors are expected to have main sequence masses in the range 
from 1 to about 8 Mq mass (the upper mass limit is not well established). As these stars 
leave the main sequence they evolve into the RGB, Horizontal Branch, AGB, and PN phases 
in the Hertzsprung-Russell (HR) diagram, before descending the white dwarf cooling track. 
The major structural changes in the post-main sequence evolution of low- and intermediate- 
mass stars occur during the RGB and AGB phases. During the RGB and AGB the stellar 
effective temperature is always lower than its main sequence value and therefore it has no 
influence on the planet's survival. However, it is during the late AGB evolution, the so-called 
thermal-pulsing AGB phase, that a planet's orbit will be most influenced, since during this 
phase the star loses most of its initial mass and reaches its maximum radius. 

The planet will spiral-in and evaporate totally (or in rare cases will accrete mass and 
become a close, low-mass companion to the star) if the planet's orbital distance is within the 
reach of the star's radius during the AGB phase (ILivio fc Sokerlll984j ). An estimate of the 
maximum planet mass that can be evaporated inside an AGB envelope can be obtained by 
equating the location of the evaporation region (where the local sound speed in the stellar 
envelope matches t he escape ve l ocity from the planet's surface) to the energy required to 
expel the envelope (ISokerl Il996l . Il998l : iNelemans &: Taurid Il998l ). The value of this maxi- 
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mum mass is very uncertain because it depends on se veral factors, such as the efficiency of 
envelope ejection (jPortegies Zwart fc Yungelson 119981) . which are largely unknown. Using 
the simplified formalism of iNelemans &: Taurid (119981 ) we find that planets with masses less 
than 0.014 Mq or 15 Mj (where Mj is the Jupiter mass) will evaporate inside the envelope 
of an AGB star with main sequence mass of 1 Mq. This mass limit is much higher, ~ 120 
Mj(0.11 Mq) well into the stellar regime, if the planet (or brown dwarf) is engulfed inside 
the AGB envelope of a 5 Mq star. 

We can therefore safely assume that planets less massive than 15 Mj that are engulfed, 
will be dissipated, because this hmit corresponds to the temperature (and therefore planet's 
escape velocity) reached inside the AGB envelope for the lowest mass stars (in the range 
examined here). Note that thi s planet mass l i mit is well below the mass of the brown dwarf 
(52 Mj) recently discovered by lMaxted et al.l (120061 ) orbiting an under-massive white dwarf. 
On the other hand, the case of total planet evaporation (for a planet in a close orbit around 
a solar- like star) has been propos ed to explain the formation of single under-massive white 
dwarfs (INelemans fc Taurislll998l ). 



The structural changes that an AGB star undergoes in response to the dissipation of 
a plane t in its i nterior are coniplex and h ave been extensively explored by ISiess &: Livio 
(Il999al lb[) and by IStruck et al.l (120021 . |200J) . The details of the destruction of such a planet 
within the stellar envelope of an AGB star are beyond the scope of the present work. 



2.1. The Evolution of the Planet's Orbit 

In the following we concentrate on Jupiter-like planets in orbits that avoid engulfment 
when the radius of the AGB star expands during the thermal pulses, and we explore the 
conditions under which the planet will remain outside the stellar surface. 

For orbital distances larger than the stellar radius, R^,, there is still a range of orbits for 
which drag effects (gravitational and tidal) are important in decreasing the planet's orbit, 
ultimately causing the planet to spiral into the AGB envelope. While the gravitational drag, 
caused by the increase in the planet mass th rough the accretion of ma terial from the stellar 



wind, is negligible for Jupiter-like planets (IDuncan fc Lissauerl Il998l ). the decrease in the 



planet's angular momentum caused by the tidal drag cannot be neglected for orbits close to 
the stellar radius during the AGB phase. 

If the planet's orbital distance is larger than the stellar radius reached when th e star 



expan ds during the AGB phase, the variation of the orbit can be approximated by (e.g. lZahn 



19771 ): 
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where r is the orbital distance and M^, is the stellar mass on the main sequence. The first 
term on the right hand side of Eq. 1 represents the increase in the orbital radius due to mass- 
loss from the star and the sec ond term rep resents the decrease due to tidal interaction, which 
is proportional to [i?*(t)/r]^ ( jZahnl 119771 ) and becomes negligible at large orbital distances. 
At small orbital distances, however, the decrease in the planet's angular mor nentum caused 
by the tidal interactiori can lead to a decrease i n the orbit by ~ 100 to 300 (ILivio fc Soker 
19831 : iRasio et al.lll996l : iRybicki fc Denisll200ll ). depending on the planet's mass and the ratio 
of the rotational to orbital angular velocity. Note however that the stellar radius during the 
AGB increases gradually. Variations in the surface luminosity and radius of an AGB star 
arise as a consequence of the thermal and structural readjustments produced during the 
after-pulse phase. In particular, i?* is not maintained after a thermal-pulse. Rather, the 
star eventually contracts as the temperature in the He-burning shell decreases and a new 
thermal pulse begins. 

The maximum stellar radius, K^"-^, is only reached briefly at the end of the AGB 
evolution as the star gets to the AGB tip luminosity. Consequently, it is safe to neglect the 
tidal term in Eq. 1 for orbital radii larger than R^°'^ . The AGB and post- AGB evolution 
of the stellar mass-loss rate, effective tempe rature, and luminosi t y that we use have been 
obtained from the stellar evolution models of IVassiliadis fc Wood! (Il993l . Il994l ). From these, 
we have calculated the maximum radii , -R™"^, reached by stars of different masses during 
the AGB phase. These values of BJ^""^ are listed in Table 1. Neglecting tidal interaction for 
r > EJ^"-^ and integrating Eq. 1, the orbital radius is then given by 



r(t) = To 



M,{t) 



(2) 



where Tq is the initial orbital radius and M^,(t) is time-dependent stellar mass caused by 
mass-loss. 

In Fig. 1 we show the evolution of the orbital radius during the late AGB phase by using 
= RT"'^ foi' each initial mass considered. In the calculation of M^,{t) we have considered 
the fact that the AGB mass-loss accounts for most of the mass lost by low-and intermediate- 
mass stars and that mass-loss previous to the AGB is negligible. Th is is a consequ ence of the 
fact that mass-loss rates on the RGB are proportional to R^L^,/M^ ( lReimerslll975l ). resulting 
in a significant amount of envelope mass lost before the AGB phase only for stars with initial 
masses M^, < 1 M©. 



The higher the initial mass, the higher the amount of mass lost during the AGB phase 
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and hence the larger the planet's orbital expansion. Planets reach a final orbital distance at 
the end of the AGB phase, determined by multiplying the initial orbit by M^^jM^D where 
Mvi/D is the white dwarf mass, given that mass-loss is negligible once the star leaves the 
AGB. The orbital expansion factors for each stellar mass are given in column (4) of Table 1. 

Another effect to consider is whether the planet becomes unbound due to the change 
in mass of the central star. Unbinding can be expected if the stellar mass-loss timescale, 
'^mass—loss-) satisfics Tjyiagg_iQgg < T/j^yni whcrc T^Q_gg_iQgg is givcu by 

Tmass—loss ~ ) 

and the dynamical timescale, r^j^n, by 

^3 

with being the stellar mass-loss rate, Mp the planet's mass, and G the gravitational 
constant. Given that Tdyn ~50yr while the shortest mass-loss timescale is Tmass-ioss ~ 10^ yr, 
it is very unlikely that a planet will become unbound due to the decrease in the stellar mass 
during the AGB phase. 

A gas planet with Mp < 15 Mj orbiting an AGB star (of any mass) at an initial orbital 
distance Vo < W^""^ will most likely evaporate inside the stellar envelope. The deuterium 
burning minimum mass limit ( 12 Mj) is traditiona lly used as the mass boundary between 



planets and brown dwarfs (e.g. ISaumon et al.lll996l ). Therefore we can safely state that gas 



planets do not survive the AGB evolution if their initial orbits lie within the reach of the 
stellar radius during the AGB evolution, whereas brown dwarfs can survive. For Tq > R™"^ 
the orbit will generally expand due to the heavy stellar mass-loss rates experienced during 
the AGB evolution. Larger differences between the initial and final mass of the star are 
experienced for the more massive progenitors, causing the orbits of planets orbiting the 
more massive stars (note that we are always refering to stars in the 1-5 mass range for 
which complete stellar evolution models exist) to be modified by the larger factors (up to 
5.5 times larger than the initial orbit). 



3. THE CHANGE IN THE PLANET'S ENVIRONMENT DURING THE 

PN PHASE 



During the post-main sequence evolution of the star the stellar effective temperature 
always remains lower than its main sequence value. However, once the star leaves the AGB 
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phase, the high mass-loss rate ceases and the remnant core moves in the HR diagram at 
constant luminosity toward higher effective temperatures into the PN stage, before the star 
reaches the white dwarf cooling track. The planet's orbit is not expected to change further 
at this stage. However, the main processes responsible for shaping PNe (high velocity winds) 
and powering PNe line emission (high stellar effective temperatures) need to be considered 
to establish the survival of a planet during this phase. 

The luminosity, mass, and ti mescale of evolution of the star during the PN phase depend 
most l y on the st ellar core mass (jPaczyhskil Il97ll : llben fc Renzinil Il983l : IVassiliadis fc Wood 
1994 : llbeni Il995l ). The stellar luminosity during this phase is within the range 3.5 to 23 
xlO^ Lq (for the lowest 0.56 and highest 0.9 Mq mass remnant, respectively) and the stellar 
temperature can reach 100 000-380 000 K (for the same core ma sses, respectiyely). The 



hydrogen ionizing photon flux, which is of the order of 10 s ( 1 Villaver et al.l l2002bl ) . is 



responsible for the PNe ionized line emission. PNe central stars also emit very high velocity 
winds (with speeds of a few thousands of kms~^) which are driven by the t ransfer of photon 



momentum to the gas, through absorption by strong resonance lines (e.g. iPauldrach et al. 
19881 ). The PN is largely shaped by the int eraction of this h i gh velocity wind w i th the 



slowly ejected material during the AGB phase (IKwok et al.lll978l : iKahn fc Westlll985l : iBalick 
19871 : iMellema fc Franklll995l : Ivillaver et aDl2002al jbh. The survival of a gas planet as the 
star evolves into the PN phase strongly depends on the planet's surface temperature, which 
ultimately determines whether or not high evaporation rates are set at the planet's surface. 



3.1. The Planet's Equilibrium Temperature 

The radiative equilibrium temperature of a planet can be estimated by balancing the 
flux received by the star to the blackbody flux re-radiated by the planet. 



^^i-Ti^J'-^) (5) 

where a is Steffan-Boltzmann's constant, Tp is the planet's equilibrium temperature, r is the 
orbital radius, is the stellar luminosity a nd A is the Bond alb edo. We have assumed an 



approximate value of A =0.5. As shown in iMarley et al.l (119991 ) Bond albedos show little 



sensitivity to the planet's gravity or effective temperature for a given stellar spectrum; the 
largest variations arise from the stellar spectral type and the presence of clouds. 

Fig. [2] shows the planet's equilibrium temperature versus orbital distance as given by 
Eq. [5] for the two extreme cases of primaries: the lowest and the highest stellar masses 
considered. Internal energy sources have not been considered in Eq. |5] since they are expected 
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to be negligible. 



A planet orbiting a PN central star will have an age in the range 0.1 to 10 Gyr as the 



star e nters the PN phase, a s sumin g that the star and the planet are coeval. iBurrows et al. 



(I2OOII ) and iHubbard et al.l (120021 ) show the predicted range of effective temperatures for 
isolated giant planets and brown dwarfs of various masses to be between 100 K and 400 K 
for the age range of interest here. The cooling theory applies equally to the structure and 
evolution of brown dwarfs and giant planets, given their convective nature. From Fig. [2] we 
see that the planet's equilibrium temperature in the presence of the PN core is much higher 
than the expected value in isolation, at orbital distances r < 20 AU and r < 50 AU for 
the 0.56 and 0.9 M© stars respectively. 



3.2. The Planet's Exospheric Temperature 

Given the high effective temperature of the star, most of the stellar flux is emitted 
at short wavelengths, for which the material in the planet's atmosphere has a high opti- 
cal depth. The incident flux will therefore be absorbed by a thin surface layer, the depth 
of which will be determined by the transparency of the material to the incident photons. 
Inward to this surface layer, the heat transfer to the interior is set by a combination of 
conduction and convection. The temperature on the planet's atmosphere has to be deter- 
mined through the energy equation considering radiative, conductive and convective terms. 
Convective effects are expected to be important in the atmosphere of a tidally synchronized 
planet, since a strong temperature gradient will be set between the day and night side of 
the planet. Howeve r, only very short period planets are expected to be tidally synchronized 



(IGuillot et al.l 119961 ) and planets orbiting PN central stars are not expected to be found at 
such short orbital periods (see §2). Therefore, it is reasonable to assume that, for the case 
under consideration, convection can be neglected as first approximation. Moreover, it has 
been shown that in the upper atmosphere of a planet (the exosphere) the density becomes 
very low and the temperature is determined mainly by the absorption of extreme ultraviolet 



and X-ray (XUV) radiation (IBauer fc Lammed l2004j ). Therefore by neglecting convective 



effects the energy equation takes the simplified form: 



V-(irVT-g.,.) = (6) 

where ^xuv is the effective radiative heat flux of the star at wavelengths shorter than 3000A for 
which the opacity is very high, T is the temperature, and K is the thermal conductivity. 
The exospheric temperature, Tf,^, can be obtained by integrating Eq. [61 This gives (see e.g. 
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Bauei]ll97ll . ll973h 



t::' - ^ p^n. (7) 

Komngaa 

where s i s the temperature depen dence of the thermal conductivity, s = 0.7 for atomic 



hydrogen (IDalgarno fc Smith! 119621 ). is the temperature at the base of the thermosphere, 
a is a factor that takes into account the rotation of the planet and varies between 0.5 and 
0.25 for rapid or slow rotators, kb is Boltzman's constant, ctc is the gas-kinetic collision cross- 
section, Ko is the thermal conductivity proportionality factor, mn the mass of the hydrogen 
atom, and g is the gravitational acceleration. The intensity of the stellar radiation reaching 
the planet's atmosphere /qo is related to gxuv by, 

gxuv = rijOaelooe'"' (8) 

where e is the heating efficiency (the fraction of the stellar flux absorbed that is transformed 
into thermal energy in the atmosphere), and r is the optical depth, with t = J rijaadz where 
rij is the number density of the absorbing constituent (atomic hydrogen), cxa its absorption 
cross section to XUV radiation and z is the atmospheric altitude. 

Since Tf, is roughly equal to the planet's equilibrium temperature, then Ti, -C T^x and the 
second term on the left hand side of Eq. [7] can be neglected. Note that in order to estimate 
Tex in Eq- cooling in the infrared has not been considered. This is justified by the fact 
that most of the molecular hydrogen (which is responsible for the cooling in the infrared) is 
expected to be dissociated in the upper atmosphere. 

In Fig. [3] we show the exospheric temperature Tex reached by a Jupiter-like planet of 1 
Mj mass due to the XUV generated by a 0.56 Mq (1 Mq main sequence mass) central star 
versus the orbital distance. The different curves are the result of changing the XUV flux 
intensity due to the evolution of the star in the HR diagram. The stellar luminosity emitte d 



below 3 OOOA has been estimated using the GMFGEN code of iHillier fc Milled fll998l . Il999r ). 
Note that for a PN central star with an effective temperature of 115 000 K (and log g = 5.7), 
97 % of the total stellar luminosity is emitted below 3 000 A. For cooler central stars (e.g. 
35 000 K), 93% of the stellar luminosity is below 3 000 A. The solid curve represents Tex 
obtained by using the stellar flux early in the evolution of the star into the PN phase when 
the stellar luminosity is = 3 600 Lq and Tefj =36 000 K; the dashed curve corresponds to 
an intermediate stage (T^, = 1 400 L0 and Tcfj =115 000 K); and the dotted curve corresponds 
to a point 10^ yr later in the evolution, when the luminosity has dropped to 100 L0 at a 
temperature of Tog =96 000 K. In the following, we explore the conditions under which an 
evaporation flow is established at the planet's surface. 
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PLANET EVAPORATION DUE TO THERMAL HEATING DURING 

THE PN PHASE 



The standard form ulation for th e gas- kinetic or thermal evaporation from celestial bodies 
was first established by I Jeand (Il925l ). Jeans's formula gives a lower limit to the rate of escape 
(by thermal evaporation) of gas from a planetary atmosphere: 



20F 



rir 



:i + E)e~ 



(9) 



where Vo = ^J2k^|T /niH is the velocity of the escaping particles, Uc is the number density of 
the escaping constituent (we assume hydrogen from now onwards) at a given critical level, T 
is the temperature, and E is the escape parameter E[R) = GMuih / RkbT where M is the 
planet's mass contained within a spherical surface of radius i?, and ki, and G are Boltzman's 
and the gravitational constants respectively. Jeans escape then usually refers to the density 
of the critical level at the ex obase, that delineates the lower limit of the exosphere (e. 



g. iChamberlain &: Hunterull987l ). The temperature of this layer, the exospheric temperature 
Tex, is much higher than the planet's effective temperature (given by Eq. [5]), and higher 
than the gas photoionization t emperature (10 000 K) and is the characteristic planetary 



temperature for thermal escape (jMoutou et al.ll200ll : iLammer et al.ll2003l : iBaraffe et al.ll2004 



Lecavelier des Etangs et al.ll2004l ). A common misconception when considering the evolution 
of a planet inside a PN, is that the evaporation rates are determined by the Jeans escape, 
with a temperature of 10 000 K fixed by photoionization. Using the ionization temperature 
the sound speed of the ionized gas in the planet's atmosphere is smaller than the typical 
escape velocity for the material and therefore no evaporation is expected. 

Equation [H] seems to imply that high evaporation rates may be achieved for high values 
of Uc- However, Jeans escape is only valid when the escape rate is relatively small and 
therefore the static structure of the atmosphere is not perturbed, that is, when the escape 
parameter E > 20-30. Note that E represents the ratio of the gravitational potential energy 
to the particle kinetic energy or, in other words, the ratio of the square of the escape to the 
thermal velocity in the upper atmosphere. For large escape fluxes, E <20, the structure 
of the upper atmosphere is modified, the gas that has escaped has to be replaced by gas 
from the lower levels thereby altering the atmospheric structure. In this case, the escape 
rate given by Eq. [HI and the thermal balance are no longer valid ( IChamberlain &: Hunten 
1987 : Bauer &: Lammer 12004). Fo r E < Ec, where Ec ~ 1.5 — 3, the atmosphere reaches 



blow-off conditions (e.g Opilj 



19631 ). The exospheric base is lifted outwards, the planet radius 



increases and the atmosphere is lost if the heat source continues to operate. Under blow- 
off conditions the planet's atmosphere will expand radially outward at relatively modest. 
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subsonic velocities, close to the surface, and will gradually accelerate to supersonic velocities 
as it moves further away. Fig. [3] shows that blow-off conditions (represented by the region 
above the solid or the dotted horizontal line depending on the Ec value assumed) are easily 
reached at the planetary atmosphere up to distances of 125 AU. 

As the stellar mass is increased, higher luminosities and effective temperatures than 
the ones used in Fig. [3] are reached during the PN phase. More intense XUV radiation 
fields will reach a planet's atmosphere orbiting a higher mass progenitor. The results shown 
in Fig. [3] therefore represent a lower limit (in stellar mass) to the development of blow-off 
conditions. Although we have assumed a hydrogen composition to compute Tex, if other, 
heavier constituents exist they will be carried away through drag forces (by the light gas 
moving at a sonic speed) once the atmosphere reaches the blow-off condition. 

The escape parameter, as well as the exospheric temperature, depend on the planet's 
surface gravity. The higher the planet's mass, the harder it is to reach the blow-off conditions 
because the particles need a much higher thermal energy to reach escape velocities at the 
planet's surface. The gray area in Fig. H] shows the orbital distance region (to 125 AU) at 
which a planet's or a brown dwarf's atmosphere will reach blow-off conditions as a function 
of the pl anet (or brown dwar f) mass. In the calculation we have used mass-radius relations 
given by lBaraffe et al. J2003h for non-irradiated extrasolar planets at 10 Gyr, and assumed 



that the exospheric temperature is given by Eq. [3 The left and right panels have been 
computed for a 0.56 Mq (1 Mq main sequence) and 0.9 Mq {5Mq main sequence) star 
respectively, early in the PN phase. We have chosen the same effective temperature of the 
star, 36 000 K, which corresponds to = 3 600 Lq and L^, = 23 000 Lq respectively. Higher 
mass stars produce more intense XUV radiation fields and blow off conditions are set at 
larger orbital distances for a given planet mass. The evolution of the star towards higher 
effective temperatures has similar effects. 



4.1. The Planet's Evaporation Rates 



Given the high exospheric temperatures reached at the planet's upper atmosphere it 
is expected that an outflow will develop as a consequence of the absorption of the XUV 
radiation. Since high temperatures can cause the outer layers of the planet to escape rapidly, 
in order to examine the planet survival, we have to estimate the evaporation rate. Note that 
at the orbital distances we are considering, t he planet will be well within the typical inner 
radius of the nebular shell (0.01 to 0.1 pc) (j Villaver et al.l l2002bl ) . so we do not expect a 
decrease in the photon flux arriving at the planet's surface due to absorption by the nebula. 
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The problem of a general outflow from a s tellar (or planetary) body can be described 
with the same set of equations used by iParkerl (119631 ) to describe the solar wind. Although 
a complete treatment of the evaporative wind requires the integration of the energy, mass, 
and momentum tr ansfer equations, we can estimate the outflowing particle flux (e.g. 
Watson et al.lll98ll ) by equating the energy input (eLxuv/4) x (i?i/r)^ to the energy required 
for hydrogen to escape GMpmn/ Rp, giving 



H 



Ar^GMpmn 



(10) 



where Ri is the planet's radius where most of the XUV radiation is absorbed, deflned as the 
level where the optical depth is unity, r is the orbital distance, Rp and Mp are the planet's 
radius and mass respectively and eL^m the fraction of the stellar XUV luminosity (see §3.2) 
that is converted into therm al energy. A proper determination of e requires the solution 
of the atmospheric structure (jYeUg 1200J) . Most calculat i ons iii the li teratu r e assu med the 



value of e = 0.63 determined for Jupiter b v IWaite et al.l (119831 ) (e.g. lYelld (120041 )). Note, 
however, that ICecchi-Pestellini et al.l (120061 ) have recently estimated the X-ray contribution 
to exospheric heating in hydrogen-ri ch planetary atmospheres and found heating rates sig- 
nificantly higher than those used by lYelld (120041 ) for planets at the same distance from the 
central star. Given the high effective temperatures reached by PN central stars we consider 
it reasonable to assume a value of e = 1. 

The determination of or ^ = Ri/Rp, is not an easy task, as it requires the full solu- 
tion of the hydrody namical escape prob lem with radiation transfer in a strongly externally 
heated atmosphere. IWatson et al.l (Il98l[) assumed a sing le layer atmosphere and an iterative 



method to solve the problem. iLammer et al.l (l2003l ) and iBaraffe et al.l (|2004| ) used Watson's 
solution for Ri which matched the observed expanded atmosphere of HD 209458b [Ri = 3 
Rj). Note that in this approach it is implicitly assumed that i?i lies below the sonic level. 

Lower mass-loss rates (by a factor of 20) than those estimated under hydro dynamical 
blow-off for HD209458b, have been obtained by lLecavelier des Etangs et al.l (120041 ). by using 
what they refer to as geometrical blow-off: a combination of Jeans escape rates and tidal 
forces. Under geometrical blow-off the high exospheric temperatures push the position of 
the exosphere to the Roche-lobe limit before hydrodynamical blow-off (or energy-limited 
mass-loss) conditions develop, resulting i n lowe r mass-loss rates than those estimated un- 
der hydrodynamical blow-off. iTian et al.l (120051 ) found an escape rate 16 ti mes smaller (for 
a heat ing e fficiency of 0.6) than the maximum evaporation rate derived by iLammer et al. 



(|2003r) and iBaraffe et all (l2004r) f or HD 209458b and consiste i it wit h that determined by 
Lecavelier des Etangs et al. (j2004 ). More recently, Jaritz et al. (2005) have shown that exo- 
planets at small orbital distances will be subjected to either hydrodynamical or geometrical 



- 12 - 



blow-off conditions depending on where the position of the exobase reaches the Roche lobe 
limit. 



In the lTian et al.l (120051 ) simulation the hydrodynamical escape of hydrogen from a plan- 
etary atmosphere was studied by allowing transonic solutions and using a two dimensional 
energy deposition layer. With this approach they found the escape rates to depend on the 
position of the heating layer as well as on the amount of energy input, with the escape rate 
exponentially increasing with the heating rate. iTian et al.l (120051 ) also found Watson's escape 
fluxes to be too high when low density hydrogen atmospheres are ir ivolved (as i n the case 



HD 209458b). At hydrogen number densities higher than 10 cm . iTian et al.l (120051 ) find 



higher scape rates than Watson. Generally, the simulations produce higher escape rates than 
the ones given by Eq. (TD] for ,^ = 3 as the density increases. This is the result of the higher 
total amount of energy absorbed in an extended atmosphere, as opposed to that absorbed 
in a single layer (the approximation used to obtain Eq. [T0|) . 

The mass-loss rates (obtained by using ^ = 3 in Eq. [TOl) from a 1 Mj planet for central 
stars (at 36 000 K) with different masses are plotted in Fig. versus the orb ital distance. 



Note that the value of ^ used in Fig. O is always smaller than that obtained by lErkaev et al. 



( 120061 ) for the exoplanets in which the mass, radius and stellar parameters are known. The 
different lines account for the different central star masses considered, with the mass-loss 
rates increasing with the stellar mass. 



It is important to note that for the geometrical hlow-off ( iLecavelier des Etangs et al. 



20041 ) to operate, the planet has to be in a very close orbit to the parent star which is 
never the case for planets to be found around PN central stars. Moreover, none of the 
calculations in the literature includes heating rates as high a s the ones expected for a plane t 
orbiting a PN central star, and as it has been shown by ICecchi-Pestellini et al.l ( l2006l ). 
the inclusion of X-ray irradiation from the star strongly increases the heating in planetary 
exospheres. Therefore, the mass-loss rates given in Fig. [5] should be considered merely as 
order of magnitude approximations to the actual mass-loss rates from a planet exposed to 
a PN central star. In addition, it is very likely that the planet will inflate as radiation is 
transformed into heat inside its atmosphere, which will lead to further increase in the planet's 
evaporation rate with our approach. An appropriate determination of the escape rate will 
require a solution of the hydrodynamic escape equations for the case under consideration. 
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5. THE PLANET EVOLUTION UNDER MASS-LOSS 

Central stars of PN do not maintain high luminosities for an extended period of time 
and therefore the evaporation rates shown in Fig. [5] will change as the star evolves. In Fig. [6] 
we show the evolution of the evaporation rate of a IMj planet orbiting a 0.56 Mq PN central 
star at two orbital distances: 1.5 AU and 5AU (solid and dashed lines respectively). The 
evaporation rate ha s been computed using £ = 3 in Eq. [TO]and the evolution of the stellar flux 



using the iPrada Moroni fc Stranierd (120021 ) models. 



from the models of IVassiliadis fc Wood! (11994r) extended to the white dwarf cooling tracks 



It is important to note that the evolution of the evaporation rate shown in Fig. E]is clearly 
an oversimplication of the real case, which involves factors other than the evolution of the 
stellar luminosity. For instance, feedback processes are important because the evaporation 
rate depends on the planet's structure which is itself determined by the planet's reaction 
to evaporation and heating. In order to establish the importance of the feedback processes 
in determining the evaporation rates we need an estimate of the relevant timescales. The 
timescale over which significant reaction takes place is given by the thermal timescale rth- 
For the case under consideration the heating by radiation is taking place in a thin layer near 
the surface of the planet where the opacity is very high. Therefore, the relevant thermal 
timescale is the one associated to the heating of the planet's external layer by an external 
heating source, r^h = GMpAM/2RL^,, where AM is the mass that is being heated and L^.. 
is the stellar luminosity arriving at the planet's surface (multiplied by a correction factor 
that accounts for the heating efficiency) . Since Jupiter-like planets are believed to be mostly 



convective, with a radiative external zone that contains ~ 0.03% of the mass (iGuillot et al. 



19961 ). it is reasonable to assume that the stellar energy is mostly deposited in this layer, 
which then gives rth ~ 10^ yr (for a 0.56 Mq star). We should now compare rth with the 
mass-loss timescale t^^, Tj^ = Mp/Mp, and if we assume that Mp = $h as given by Eq. [TO] 
then ^ 10^ yr (for 1 Mj and 0.56 M©). Therefore, for range of orbital distances, Tj^-j- < rth 
which implies that the planet's mass is decreasing faster than the time the planet has to 
readjust to a new thermal equilibrium. 

The planet has two ways to adjust to a new thermal equilibrium: to increase its radius 
or to increase its effective temperature. If the planet reacts to the external heating by 
increasing its effective temperature, then the evolution of the mass-loss rates computed 
using Eq. [TU] and the time- dependent stellar flux are a good first order approximation to the 
actual evaporation rates. However, given the high mass-loss rates involved, the planet may 
react by increasing its radius on a timescale given by rth, in which case the mass-loss rates 
will be larger than estimated. In fact, once the mass- loss timescale becomes much shorter 
than the readjustment (thermal) timescale (r^^-^ < O.lrth), a runaway mass-loss may ensue 
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e.g. ISaraffe et alJl2003h 



5.1. Other Evaporative Processes 



There are other non-thermal processes that might contribute to the planet evaporation, 
e.g. photodisociation and ram pressure stripping. An increase in the planet's radius will 
decrease the density of the outer layer allowing the ionization front to propagate inwards, 
deeper into the structure of the planet. It can be easily shown that given the high densi- 
ties involved, photoioniziation in itself cannot destroy a Jupiter-like planet. However, the 
ionization of the gas increases the pressure (as it increase the number of particles and the 
temperature), causing an expansion of the outer layers that is driven from the outside and 
not from the luminosity of the planet's interior. It has even been suggested that a Jupiter 
like planet could be detected by the variable hydrogen recombination line emission it woul d 
emit in its outer atmosphere which is photoionized by the PN central star (jChu et al.ll200ll ). 



Ram pressure stripping caused by the stellar wind might become a very efficient process 
in removing mass from the outer layers, once the planet expands. Ram pressure stripping 
could dominate during the early stages of the evolution of the PN phase, when high velocity 
winds (with still relatively high densities) directly impact the planet's surface (at that stage 
PwvIj \staT > {Mp/4:TTRp)ve \pianet with the plauet's scape velocity). Ram pressure stripping 
due to the stellar pr oper motion has been shown to be a very efficient process in removing 
mass from PN shells (jVillaver et al.ll2003l : IVillaver fc Stanghellinill2005l ). Furthermore, evap- 
oration rates might increase significantly when one takes into account the rotation of the 
planet JBurkelll969h . 



As the star evolves, the wind velocity increases, reaching a velocity above which the 
shocked gas is not able to cool down radiatively. An adiabatic shock then develops at the 
interaction region between the high-velocity stellar wind and the dense material ejected pre- 
viously during the AGB, forming what is known as the hot bubble. The post-shocked gas 
temperature can easily reach 10^- 10'' K and h as been observat i onaly confirmed through th e 
detection of diffuse extended X-ray emission (IChu et al.ll200ll : iGuerrero et al.l |2002| . |2005| ). 
The outer radius of the hot bubble c an reach "^0.1 pc a nd the reverse shock region is ex- 



pected to be very close to the star (jVillaver et al.ll2002bl ). Therefore, the planet's orbit will 
be immersed inside an environment at the hot bubble temperature of 10^-10'' K. The hot 
bubble, however, does not develop inmediately (jVillaver et al.ll2002bl ). because the wind ve- 
locity needs to be high enough to develop an adiabatic shock. It is then expected that the 
temperature of the planet's atmosphere will be first set by the absorption of XUV radiation 
as described in §3.2, and only during later stages of the star's evolution the hot bubble will 
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provide an additional heating source. 

It can be expected that photoionization, ram pressure stripping and the presence of the 
hot bubble will increase the evaporation rate from a gas planet orbiting a PN central star. 
All these processes have not been included in the evaporation rates given in Fig. [5] (for ^ = 3 
which should therefore be regarded as representing conservative lower limits. 



6. PLANET SURVIVAL 

We have shown that for small orbital distances the planet's structure will change on 
timescales longer than the mass-loss timescales and that the most likely response will be an 
increase in the planet's radius. As a first approximation and since we are not computing 
the rate of change of the planet's radius (a problem that requires a full hydrodynamic 
calculation) we have estimated the planet's survival by integrating the evaporation rates 
obtained by evolving the stellar flux under different assumptions for ^. As mentioned in §4.1, 
^=3 represents the most conservative scenario, most likely a lower limit to the evaporation 
rate of the planet. Note that for this particular problem, an increase in the planet's radius 
will naturally lead to lower densities in the outer layer, making the planet more vulnerable 
to the erosion caused by photoionization and ram pressure stripping. Since we have not 
quantified all these effects when estimating the evaporation r ates in Eq. [TDl using .^=10 is 



probably a reasonable assumption. Note that lTian et al.l (120051 ) have shown that the position 



of the exobase for the extended atmosphere of HD 209458b is close to ^ = 10. 

Table 2 lists the percentage of the planet that is evaporated until a 0.56 Mq star enters 
the white dwarf cooling track. Column (1) gives the initial planet mass, column (2) the orbital 
radius, and columns (3), (4) and (5) give the percentage of the planet mass evaporated by 
using .^=3, 5, and 10 respectively. Note that although we have extended the evolutionary 
time into the white dwarf cooling sequence, most of the planet's evaporation takes place as 
the star is evolving during the PN phase. 

In Fig. [7] we show the region (on the planet mass versus orbital distance plane) inside 
which Jupiter-like planets will be destroyed, as the star evolves off the main sequence. The 
left and right panels are for planets orbiting a 0.56 and 0.9 Mq star respectively (which 
correspond to 1 and 5 Mq main sequence masses). The light gray area represents the region 
inside which the planet will be engulfed and most hkely destroyed as the star ascends the 
AGB. The dark gray and shaded areas (computed for ^=3 and 10 respectively) represent the 
regions for which planets will lose 50 % of their mass before the star enters the white dwarf 
cooling track due to the evaporation caused by thermal heating. Note that the region for 
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which the planet will undergo total evaporation due to thermal heating is inside the AGB 
stellar radius for the 0.9 Mq star. 

Giant planets are thought to be formed in cool regions of the protoplanetary disk- 
beyond t he ice line- where there is enough solid mat erial to produce a core th at will capture 



the gas (ISasselov fc 



19981 : lArmitage et al 



iCcar 



Eooi) i 



2000 : iKornet et al.ll2004l ). Planet migration (e.g. iTrilling et al. 



is then used in order to explain the close dis t ances at which such 



planets are observed aro und sun-like stars (e.g. iMayor fc Quelog|l995l : IVidal-Madjar et al. 
20041 : iMarcv et allboosi ). Most of the planets recently discovered are giant planets orbiting 
sun- like stars at orbits < 5 AU. 

The conditions for planet survival as the star evolves off the main sequence depend on 
the initial mass of the star. From Fig. 7 we see that most of the close-in planets (r < 1.5 
AU) will be destroyed as they get engulfed by the star during the AGB phase (see Table 1 
for a list of the maximum stellar radius reached during the AGB versus the stellar mass). 
As the star leaves the AGB and enters the PN phase, high effective temperatures at very 
high luminosities set up an evaporation flux at the surface of the planet. At certain orbital 
distances the evaporation rates are high enough to cause a total destruction of the planet. 
By integrating the evaporation rates as the star evolves during the PN into the white dwarf 
cooling track, we find that Jupiter-like planets will be destroyed if they remain at orbital 
distances r< 5 AU from a low mass white dwarf (Mwd < 0.63 Mq) and large planet ablation 
is expected up to 10 AU. In particular, Jupiter in our own Solar system is barely expected 
to survive. More massive stars evolve very fast during the PN phase and do not maintain 
high evaporation rates long enough to cause planet destruction, unless the planets were to be 
found at small orbital distances (r<2.5 AU). However, we have shown that planets orbiting 
the more massive PN central stars cannot be found at small orbital distances: if a planet 
orbiting a 0.9 MQprogenitor is to survive AGB engulfment then its orbit has to be at r> 29 
AU. 

Over the last few years several groups have surveyed white dwarfs in the search for 
planets with no posi tive results reported so far. The wide field proper motion survey around 
261 white dwarfs by lFarihi et al.l (120051 ) was sensitive to 100 to 5000 AU with masses down 
to 52 Mj, and their deeper near fie ld search of 8 6 whit e dwarfs could detect 10 Mj at 
separations between 50 and 1100 AU. lMuUally et al.l (120061 ) have recently surveyed 124 white 
dwarfs with Spitzer an d found no evidence for planets with M > 5 Mj. Using direct imaging 
Zinnecker et al.l (120061 ) have ruled out the presence of giant planets (6-12 Mj) around the 7 
known white dwarfs of the Hyades cluster. A plan et (4 Mj) was discovered in the binary 
system Gliese 86 by iMugrauer &: Neuhauserl (120051 ) orbiting the Kl dwarf in Gliese 86A, 
which has a 0.55 M0mass white dwarf companion, Gliese 86B, at 21 AU. At this distance 



-17- 



from the white dwarf, we expect the planet to have been losing mass during the PN phase 
at such a low rate that it would not have affected its survival. 

On the other hand, three white dwarfs are known to posses a n infrared excess, GD 362 



fIZuckerman k BecklinI 119871 : iBeckhn et all 20051 : iKilic et al.ll2005r ). G 29-39 (Jura 



20031 ) and 



the recently discovered GD56 (IKilic et al.l 120061 ) . To explain G 29-38 lJural (120031 ) proposed 
a disk formed by the tidal destruction of an asteroid. The destruction and accretion of a 
planetary body continues to be the widely accepted scenario to explain the high anomalous 
photospheric abun dances in the massive white dwarf G D 362, where a massive accretion disk 
has been detected (IBecklin et al.ll2005l : iKilic et al.ll2005l ). We have shown that planets around 
the most massive white dwarfs are expected to be found at r> 30 AU, due to the expansion 
of the planetary orbit during AGE mass-loss. How a planet at such a large distance from 
the star would then be disrupted (in order to form an accretion disk) is by no means clear. 
Another more plausible scenario for the presence of accretion disks a round mass i ve wh ite 
dwarfs, is that they are formed by the merger of two white dwarfs (ILivio et al.ll2005l ). a 
scenario that can also explain the presence of planets around massive white dwarfs at small 
orbital distances. 



7. CONCLUSIONS 

Planets at orbital distances within the reach of the stellar radius during the AGB phase 
will spiral-in and totally evaporate. The AGB stellar radius depends on the stellar mass and 
therefore so does the distance from the parent star at which a planet is destroyed, r< 1.6 
AU for M,=0.56 Moand r <5.3 AU for M,=0.9 Mq. 

If the planet avoids engulfment its orbit will typically expand. The final masses of white 
dwarfs are primarily established during the AGB phase due to heavy mass-loss. The more 
massive progenitors (when on the main sequence) lose larger amounts of mass during the 
AGB, thus the planets orbiting these stars will experience the largest orbital readjustments. 

We show that as the star evolves into the PN phase, high temperatures are reached in 
the outer atmosphere of the planet resulting in the development of blow-off conditions. As 
the star evolves through the PN phase into the white dwarf cooling track: 

1. An general outflow is established at the surface of a substellar object. The evaporation 
flux decreases with the orbital distance, and as the star evolves, with higher evaporation 
fluxes expected for planets orbiting the more massive stars. The evaporation rates are 
sufficiently high (10~^ Mjyr~^ at 5 AU for a IMj planet) that these evaporating planets 
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may be det ected using spectros c opic o bservations, similar to those performed for HD 
209458b bv lVidal-Madiar et al.l J2003h . 



2. Under these evaporation rates, planets with Mp <1 Mj will not survive the PN phase 
if located at orbital distances r < (3-5) AU. 



3. Planets with Mp > 2Mj survive the PN phase down to orbital distances of 
around low-mass (Mwd ~ 0.56 Mq) central stars. 



3AU 



Planets around white dwarfs with masses Mwd ^ 0.7 M0 (that formed from single 
stars with masses Mms ^ 2.5 M0 ), are generally expected to be found at orbital radii 
r > 15 AU due to the effects of mass-loss on the AGB phase. If planets are found at 
smaller orbital radii around such massive white dw arfs, they had to form as a result of 
the merger of two white dwarfs (ILivio et al.ll2005l ). 



Acknowledgment: We are grateful to James Herald for providing us with the stellar fluxes 
from the CMFGEN code. 
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Fig. 1. — Evolution of the orbital distance caused by the stellar mass-loss during the AGB 
phase. The time is given in logarithmic scale, with to representing the begining of the AGB 
evolution. Each curve represent a different PN central star mass, marked on each curve and 
with different line-styles. We have used an initial orbital distance that equals the maximum 
radius reached by the star during the AGB phase. 
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Fig. 2. — Left: the solid line represents the planet's equilibrium temperature (from Eq. [5]) 
versus orbital distance for a star of 0.56 Mq as it enters the PN phase. The dashed lines 
follow the evolution of the star up to 1.1x10^ yr after the PN phase starts. The lowest 
equilibrium temperature corresponding to the latest point in the evolution. Right: the same 
for a 0.9 Mq star evolved to 1 xlO^ yr into the PN phase. 
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Fig. 3. — Exospheric temperature (in Logarithmic scale) reached at the atmosphere of 
a 1 Mj planet orbiting a 0.56 Mq star (1 Mq main sequence) plotted versus the orbital 
distance. The different curves (solid, dashed, and dotted) are the result of changing the 
XUV intensity due to the evolution of the star and are chosen at luminosities of 3.6, 1.4 and 
0.1 xlO^ Lq respectively which correspond to stellar effective temperatures of 36000, 115 000 
and 96 000 K respectively. The solid horizontal line represents Tgj,, or the escape parameter 
Ec—3 (as defined in the text), above which the atmosphere reaches blow off conditions. The 
position of Ec—1.5 is also shown as a dotted hne. 
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Fig. 4. — The grey area represents the parameter space region (planet's mass versus orbital 
distance) for which the planet's atmosphere reaches blow off conditions using a critical escape 
parameter, E^. = 3. Left: computed for planets embeded in the PN environment created 
by the lowest mass PN progenitor 0.56 Mq (1 Mq main sequence mass) at an early age, 
when L*=3.6xl0^ and Tg// = 36000 K. Right: the same but for a 0.9 Mq central star. 
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Fig. 5. — Mass-loss using ^ = 3 (see text) (in logarithmic scale and Mj yr""^) versus orbital 
distance of a Jupiter-like planet under hydrodynamic escape conditions. The solid line is 
for the radiation field emited by the lowest mass core 30 000 yr after the star enters the PN 
phase of evolution. The different lines are for all the other stellar masses considered with 
the radiation (and therefore the mass-loss rates) increasing with the stellar mass. 
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Fig. 6. — The solid line represents the evolution of the evaporation rate (in logarithmic 
scale and Mj yr~^) of a 1 Mj mass planet under hydrodynamic escape conditions and at 
a distance of 1.5 AU from a . for a 1 Mostar. The same for the dashed line but for a 
distance of 5 AU. Time zero is set as the star leaves the AGB phase . The evolution of the 
st ellar luminosity during the PN phase has been taken from the stellar evolutionary models 
of IVassiliadis fc WoodI (119941) and it has been extende d to the white dwarf cooling sequence 
using the models of iPrada Moroni fc Stranierd (|2002| ). 
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Fig. 7. — Left: The dark gray area represent the region where a planet loses 50 % of its 
mass when orbiting a 0.56 M0 star by using ^ = 3 (see text). The same for the region with 
parallel line filling but using ^ = 10. The light gray area represents the orbital distance 
reached by the stellar radius during the AGB phase. Right: the same for planets orbiting 
0.9 M0 central star. 



Table 1. MASSES, ORBITAL DISTANCES 



Mms 


MwD 




r{t)/ro 


[Mo] 


[Mo] 


[AU] 




1.0 


0.56 


1.59 


1.8 


1.5 


0.60 


2.43 


2.5 


2.0 


0.63 


3.02 


3.2 


2.5 


0.68 


3.33 


3.8 


3.5 


0.75 


4.18 


4.9 


5.0 


0.90 


5.25 


5.5 
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Table 2. % EVAPORATION FOR 0.56 Mq STAR 



Mp[Mj] 


r[AU] 


e=3 




e=io 


1 


1 


100% 


100% 


100% 




2 


30% 


100% 


100% 




3 


20% 


40% 


100% 




4 


10% 


20% 


100% 




5 


5% 


20% 


70% 




10 


1% 


5% 


20% 


5 


1 


4% 


10% 


52% 




2 


1% 


3% 


10% 




3 


0.4% 


1% 


5% 




4 


0.2% 


0.6% 


3% 




5 


0.1% 


0.4% 


1.6% 




10 


0.04% 


0.1% 


0.4% 



